Quantum Computation Using Vortices and Majorana Zero Modes of a p^: + ipy 
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We propose to use the recently predicted two-dimensional 'weak-pairing' p^ + ipy superfluid state 
of fermionic cold atoms as a platform for topological quantum computation. In the core of a vortex, 
this state supports a zero-energy Majorana mode, which moves to finite energy in the corresponding 
topologically trivial 'strong-pairing' state. By braiding vortices in the 'weak-pairing' state, unitary 
quantum gates can be applied to the Hilbert space of Majorana zero-modes. For read-out of the 
topological qubits, we propose realistic schemes suitable for atomic superfluids. 
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Introduction. Topological quantum computation re- 
quires particles that have non-Abelian statistics under 
interchange and braiding. Under pairwise interchange 
of particle coordinates, the many-body wave-function of 
particles following non-Abelian statistics transforms via 
a unitary transformation in the Hilbert space of a de- 
generate set of wave- functions. Such particles can arise 
as the low-energy excitations of a topological phase of 
matter. One such system, which has been recently dis- 
cussed in this context [ll l^j 13, Ij) la] j is the experimentally 
observed 0, = | fractional quantum Hall state of 
a two-dimensional (2D) electron gas, where v is the fill- 
ing fraction of the electrons. Another promising system 
is a spin-triplet (5 — 1), 2D, + ipy superconductor, 
in which certain vortex excitations have zero energy Ma- 
jorana modesj^ in their cores, which endow these vor- 
tices with non-Abelian statistics llll . For a spin- 
polarized (spinless) Px + ipy superconductor, ordinary 
vortices with vorticity A'^ = 1 have such Majorana modes 
bound at the core 12]. For a, Px -\' ipy superconductor 
with Sz = 0, analogous to the A phase of He 3 0|, only 
the higher-energy vortices with N — ^ - not the low- 
est energy A^ = 1 vortices - have the Majorana modes. 
Nevertheless, it is possible to quench the spin-orbit en- 
ergy, which acts as a confining potential between two 
N = ^ vortices, by applying a magnetic field Q, and 
thus to make the Majorana modes potentially realizable 
in experiments. Based on this appealing idea, and the 
prospect that strontium ruthenate may be a quasi-2D, 
Sz = Q,Px-\- ipy superconductor [la,llall3) it has been 
recently proposed [iSj to use thin films of this material as 
a system which realize non-Abelian statistics of quasipar- 
ticles, associated with these unusual half-quantum (i.e., 
N = ^) vortices. 

Although there is nothing in principle to invalidate 
such a strategy, practical diSiculties may arise due to the 
lack of quantum coherent motion of vortices in the films. 



Moreover, since one needs to apply a threshold magnetic 
field to quench the spin-orbit energy, there will be a rel- 
atively high concentration of the half-quantum vortices, 
thereby rendering independent braiding experimentally 
challenging. Finally, and most importantly, since the 
quasiparticles of a superconductor are chargeless, and 
the Majorana modes are also spinless, there is no sim- 
ple way to couple to the state of a qubit after a braiding 
operation has been performed. This makes reading out 
the state of the qubit difficult. Hence, even though the 
very realization of non-Abelian statistics through the ob- 
servation of these vortices is an exciting goal in itself, 
and for the purposes of topological quantum computa- 
tion several ideas to overcome the difficulties mentioned 
above were recently proposed it will really pay to 
have a Px + ipy superfluid system where vortex motion 
is likely to be coherent, A^ = 1 vortices themselves have 
non-Abelian statistics so that their concentration can be 
independently kept low, and a natural read-out scheme 
exists. 

With the recent observation of a p-wave Feshbach res- 
onance in spin-polarized '*°K and ^Li atoms in optical 
traps 19, 20, 21], just such a system - an 'artificially' 
created Px+ipy superfluid of spinless fermions - may now 
be within experimental reach. Exotic non-Abelian statis- 
tics is thus tantalizingly close to fruition in these systems. 
Since the atoms are in identical spin states, s-wave scat- 
tering is Pauli-prohibited and a p-wave resonance domi- 
nates, allowing the tunability of the atom-atom interac- 
tion in L = 1 channel. Recently, it has been theoretically 
shown jl^, that such interactions have the potential 
to realize various p-wave superfluid states, among them, a 
Px+iPy state in the so-called 'weak-pairing' phase (chem- 
ical potential /i > 0) in both three and two dimensions. 
As a function of the Feshbach resonance detuning, which 
controls fi, this phase undergoes a topological quan tum 
transition to the strong-pairing phase (/i < 0) [22]. In 
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2D, the phase with ^ > is topologically non-trivial be- 
cause it supports zero-energy Majorana modes at vortex 
cores while, as we show below by explicitly con- 

structing the zero-mode wave- function, they disappear in 
the topologically trivial strong-pairing phase. The weak- 
pairing phase, then, is suitable for use in the hardware of 
a quantum computer. Since ordinary vortices themselves 
exhibit non-Abelian statistics in this case, they can be 
created at a low density, allowing, in principle, indepen- 
dent braiding. These vortices are also expected to be 
light due to the high degree of coherence possible in op- 
tical traps, and so it is much easier to maintain quantum 
coherence during the braiding operations. Finally, as we 
discuss later, since the atoms, unlike the electrons in a 
superconductor, have internal energy levels, this inter- 
nal structure can be manipulated to read out the states 
of the qubits after the braiding operations perform the 
quantum computation. 

The weak and strong pairing phases and the fate of 
the zero mode. The BCS Hamiltonian for a system of 
spin-polarized (spinless) fermions in a two-dimensional 
spin-triplet p-wave superfluid state is given by, 

where ^^{x) is a two-component column vector, ?/'(^) = 
(c^(x), c(x))"'", and TC{x,x') is the matrix, 

nix,x') = i^-^)Six-x')a^ + ^^a+-^^i^a- 

(2) 

Here, m is the fermion mass, fj, is the chemical potential, 
and A(a;, x') is the gap function. We take ?i = fcs = 1 in 
this paper. In a uniform p^ + ipy state, the gap function 
takes the form in momentum space, A{p) = ^{px+ipy)- 
For n < 0, the fall-off of the pair relative wave func- 
tion g{xi — X2) is exponential (pairs are tig htly bound), 
whereas, for ^ > , it is algebraic Following 
Ref. we identify the system to be in the strong- 

pairing phase for < and in the weak-pairing phase 
for /i > 0, the two phases separated by a topological 
phase transition |24| . The weak-pairing phase supports 
zero-energy Majorana fermions at the vortex cores. By 
explicitly constructing the bound state wave function us- 
ing the Bogoliubov - de Gennes (BdG) equations, below 
we show that, for /i > 0, the zero-mode quasiparticle is 
self-hermitian (Majorana), and, as /i is tuned to nega- 
tive values, these modes disappear in the topologically 
trivial strong-pairing phase. For p^ + ipy wave supercon- 
ductors, an alogous BdG equations have been discussed 
in Refs.iillllll 

To construct the eigenfunction of the Hamiltonian, 
Eq. 121 for the zero-energy state, if any, in the presence 
of a vortex, we model the vortex by assuming the gap 
function to be zero inside a circular area of radius ^ (co- 
herence length). Outside this radius, the gap function 



takes the form A{p) = exp{i9/2){px: — ipy) exp(i0/2), 
where the total order parameter phase, 9, rotates by 27r 
around the vortex with unit vorticity. In polar coordi- 
nates (p, 6), for p < the BdG equations take the form, 

1 Id 1 52 

where (j){x) — {u{x),v{x))'^ . For the zero-energy state, 
we take the angular momentum operator I = —i-§Q to 
have eigenvalue zero. The remaining parts o f Eg . El imply 
just ordinary Bessel equations of order zero |28| for both 
u and V. Since one of the two independent solutions is 
divergent at the origin, we find the solution for 0, 

0(p) = ^Jo(v^P)C, (4) 

where Jo is the Bessel function of the first kind of order 
zero, A is a constant, and C is a constant spinor. 

For solutions with p > we note that the 

gap operator can be written in polar coordinates as, 
-z^ exp(-z0/2)(^-i^) exp(z0/2) = _zA.(|. + _L)_ 

Using this, and for zero angular momentum, the 
BdG equations for the zero-energy state, on multiplica- 
tion by —2m<7z, can be written as, 

(5) 

The solutions to this equation which are well- 
behaved at /9 00 can be written as 4>{p) = 
Xip) exp(-^/3)(l, -i)^, where xip) satisfies 

[|l + i|. + (2mM-^)]x(p)=0, (6) 
op'' p op Vp 

where vp = This is again Bessel equation of order 
zero. Since both solutions are well-behaved at infinity 
(they are asymptotically sinusoidal) , the general solution 
for (pip) can therefore be written as, 

cj^ip) = [BMkp)+CY,{kp)] eM^j~—p)ih (7) 

4 VF 

where Yq is the Bessel function of the second kind of 
order zero, k — \j2mp — A'^/v'^, B, C are constants, 
and the phase factor e*~ is for equal distribution of phase 
between (j) and 0^ (see below). Next, to get a complete 
solution for the zero-energy state, one needs to match 
the wave-function and its derivative at p = ^, and also 
normalize the function in all space. These conditions will 
provide three equations for the three constants A, B, and 
C, which can then be straightforwardly solved in terms 
of the known parameters. Once the solution 4'{p) and, 
in turn, u{x),v{x) are known, the quasiparticle operator 
for the zero-energy state can be written as 

-fl^ i d^x{u{x)c\x)+v{x)c{x)) (8) 
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For the overall phase choice e'^ in Eq. Ul one can see 
from Eq. |S1 that 7o = 7o : the zero energy state is a 
self-hermitian Majorana state. 

To see what happens to the zero-energy state in the 
strong-pairing phase, this phase can be accessed, in the 
spirit of Ref. IJ, by staying within the mean- field BCS 
theory and taking /i < 0. Ref. Q argued that the edge 
of a vortex could be viewed as a wall separating vacuum 
with /i large and negative inside the core from the conden- 
sate outside. Therefore, in the weak-pairing condensate 
only, the edge acts like a domain wall between strong 
(inside the core) and weak-pairing phases, while, in the 
strong-pairing condensate, nothing interesting should re- 
sult in the core. By extending the mean-field theory to 
fi < ,2^^], and replacing fi by Eq. |31(for zero 

angular momentum) and Eq.^now imply modified Bessel 
equations of order zero for inside and outside the core, 
respectively. The solution 4>'{p) is now given by only 
one of the two modified Bessel functions of order zero in 
each case, Io{y/2m\fi\p) for p < and Kq{k'p) for p > £ , 
since the other one is divergent in the relevant region |28l | . 

Here, k' ^ {2m\p\ H — -r)^ ■ The corresponding constants 

multiplying the solutions drop out when one matches the 
solutions and their derivatives at p = ^, and divide one 
equation by the other. For generic values of the parame- 
ters, the resulting equation does not have a solution [j^ . 
and therefore, we do not expect a zero-energy state in 
the strong-pairing phase. 

For the sake of completeness, we mention here, that 
even for a spin-triplet superconductor with = 0, in 
the weak-pairing phase, there are two zero-modes at the 
core of a vortex with = 1, one for each quasiparticle 
spin. However, these modes are not Majorana modes, 
since there is mixing of up and down spin 'c'-operators 
in the definition of the quasiparticle operator analogous 
to Eq. |S1 Moreover, in the presence of any spin-flip scat- 
tering, these degenerate modes will mix and split. This 
is why one has to consider the exotic N = ^ vortices 
in these systems in an attempt to realize particles with 
non-Abelian statistics 18] . 

Non-Abelian statistics and unitary operators in the 
Hilhert space. When the system is in the weak-pairing 
superfluid phase, a dilute gas of vortices can be cre- 
ated. Suppose there are 2n such vortices in the opti- 
cal trap. Each vortex will have a zero-energy Majorana 
fermion attached to the core. For 2n vortices, there 
are 2n such fermions, which we denote by 7^, where 
i counts the vortices. The Majorana fermions can be 
combined pairwise to create n complex fermionic states, 
Q = l2i + «72i-i,cJ = 72i - i72i-i- Each one of these 
complex fermionic states can be either occupied or unoc- 
cupied, giving rise to 2"-fold degeneracy in the Hilbert 

space protected by the gap, o^o ~ 7^, where Aq is the 
amplitude of the pairing gap and e^^ is the Fermi energy, 
to the first excited state in the vortex core. Two states 



of a representative qubit are identified with the absence 
(|0)) or the presence (cj|0)) of a superfluid quasiparticle 
in the fermionic state constructed from 72^-1 and 72;. 
Note that the two states are degenerate and are not di- 
rectly associated with any particular vortices. It is this 
non-locality that protects the qubits from decoherence 
due to the environment, which acts through local opera- 
tors. For initialization of the qubits, note that creating 
vortices in pairs from the vacuum will always put each 
pair in the |0) state at zero temperature (T). At finite 
T, there is always a non-zero probability that a fermionic 
quasiparicle will end up on a vortex pair. Since we can 
read these non-destructively (see below), we can correct 
or discard the |l)'s. We now briefly describe the physics 
[loj behind the unitary transformations in this space, in- 
duced by braiding of the vortices around one another. 
These unitary transformations can be fruitfully utilized 
for quantum computation. 

Using the property that 7's carry odd charge with re- 
spect to the gauge field of a vortex with unit vorticity, 
that is, 7 —7 for a phase change of 27r, it follows 
that, upon interchange of two neighboring vortices 1 and 
2, 71 — > 72, but 72 — > —71 The unitary operator 

in the two-dimensional Hilbert space that enforces this 
transformation is given by, 

Ti = exp( J7271) = exp(z J(2ctc - 1)), (9) 

where, c = 71 -|- 1^2- This operator can be written as a 
2x2 matrix in the space of states spanned by |0) and 
1 0) . Likewise, for four vortices 71 , 72 , 73 and 74 , the uni- 
tary braiding operators can be written as 4 x 4 matrices 
in the space spanned by the basis states |0), c||0), cj|0) 
and c|c2|0), where c\ = 71 -t- 172, and c\ = 73 + *74- 
In the case of 2n vortices, the braiding operators are 
2" X 2"-dimensional matrices: they form a matrix rep- 
resentation of the braid group in two dimensions. Upon 
braiding of two vortices, an initial state, which is now a 
2"-dimensional vector in the space of degenerate states, is 
multiplied by these matrices and gets transformed to an- 
other vector in this space. It is these unitary transforma- 
tions that can be utilized to build unitary quantum gates, 
and this is the essence of topological quantum computa- 
tion. Note that T needs to be kept lower than uq. For 
Feshbach resonance-superfluids, which are near the BCS- 
BEC transition (recall that it is actually a transition be- 
tween the weak- and strong-pairing phases, unlike in the 
s-wave case, in which it is merely a crossover), Aq ef, 
and so T should be kept much lower than e^, which is 
realizable in these systems |3^. For this method to suc- 
ceed, it is imperative that the vortices can be braided 
around one another like independent particles, which re- 
quires a low density of vortices, and their movement be 
quantum coherent, both of which are achievable in opti- 
cal traps. Below, we show that the atoms in optical traps 
also offer a natural strategy for determining the state of 
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a qubit after a computation has been performed. 

Reading out the states of the qubits. A central ques- 
tion in the above scheme is how to determine the state 
of a qubit after a computation has been performed. The 
two states of the qubit, as we described above, are distin- 
guished by the presence or absence of a superfLuid quasi- 
particle at the complex fermionic state when two vortices 
are fused together. However, since these are quasiparti- 
cles of a superfluid, they are chargeless. Moreover, since 
the Majorana fermions are spinless, the complex fermion, 
which is a linear combination of two Majorana fermions, 
is also spinless. Therefore, one does not expect an excess 
of charge or spin due to the presence of a quasiparticlc 
at the core of the composite vortex. One may look for 
subtle differences in charge or energy distribution at the 
core due to the presence or absence of a quasiparticlc, 
but these may be experimentally difficult to achieve. 

A different approach, suitable for atomic (or molecu- 
lar) superfluids only, is to use the internal energy levels 
of the atoms themselves. The basic point is that, if there 
is an unpaired atom at the core of the composite vor- 
tex (the qubit is in the state c'''|0)), photons from a laser 
can be absorbed to excite the atom to an appropriately 
chosen excited level. If there is no quasiparticlc there 
(the qubit is in the state |0)), there will be no absorp- 
tion at this frequency. Note that, during this process, 
one might end up exciting Cooper pairs from outside the 
core as well. However, to excite an atom bound in a 
Cooper pair with another atom in an identical internal 
state, one first needs to break the pair, costing an energy 
2Ao. Thus, from this process, photons can only be ab- 
sorbed at a frequency shifted by 2Ao. Since the typical 
spontaneous emission rate, ~ 0(27r x 1 MHz), in such a 
detection process is much larger than Aq ^ 27r x 11 KHz 
[sof. this method can be applied only through interme- 
diate states which induce a much larger energy splitting 
between paired and unpaired atoms. Here we illustrate 
this reading out scheme using ^"if atoms, although the 
technique is applicable to other species as well. 

Suppose the atoms in the superfluid are in the 4^iS'i/2 
hyperfine ground state \i) = If = 9/2, mp = —7/2) in 
the case of p-wave resonance pj( . To determine whether 
there is an unpaired atom inside a composite vortex, 
a two-photon Raman pulse is applied that transfers 
the unpaired atom to another hyperfine state \j) = 
|F = 7/2, TTiF = — 5/2). The frequency difference be- 
tween the two Raman lasers is adjusted to be resonant 
with the hyperfine splitting between states \i) and \j) for 
the unpaired atom, but has a 2Ao detuning for paired 
atoms due to the energy cost to break a pair. The lasers 
have maximal intensities located at the core of the vortex 
and their beam waist width w ~ 1.5 fim is much smaller 
than the typical distance (> 10/im) between vortices [sij . 
allowing individual access to the qubits. The Rabi fre- 
quency of the Raman pulse is chosen to have a Gaussian 
shape 51 = fig exp (— Wpf^) (— t/ < t < tf) to reduce 



the impact on paired atoms '32'|. For a set of parame- 
ters Aq = 27r X 11 KHz, wq = Ao/2, and flo = l-77wo, 
tf ^ 5/ujo we find that the unpaired atom is com- 
pletely transferred from state \i) to \j) by the Raman 
pulse, while the probability for the paired atoms to be 
excited to state |j) is about 6 x 10^^ and may therefore 
be neglected. 

To obtain a cycling transition necessary for the de- 
tection of the unpaired atom, tt Raman pulses are ap- 
plied to transfer the unpaired atom to state \k) = 
\F — %/2,mp = 9/2). Because of large Zeeman split- 
ting between different magnetic sublevels, these Ra- 
man pulses may be performed in a short period (no 
longer than 100/xs). A focused cr "'"-polarized detection 
laser resonant with the cycling transition \k) — > |^) = 
|5^P3/2 : F — 11/2, mp = 11/2) is then applied to detect 
atoms at state \k). Here we choose the 45' — > 5P instead 
of 45* 4P transition for the detection laser to obtain 
smaller diffraction limit as well as smaller spontaneous 
decay rate. In the experiment the magnetic field 
B « 200 G for the p-wave Feshbach resonance, which 
yields an effective detuning 5 « 27r x 170 MHz for the 
paired atoms at state \i). The ratio between the num- 
ber of the spontaneously emitted photons by paired and 
unpaired atoms is estimated to be {T/25)'^ « 1.2 x 10~^, 
where F « 27r x 1.2 MHz is the decay rate for the ex- 
cited state Therefore, in the detection process, the 
impact on paired atoms may be neglected, and resonant 
fluorescence is observed if and only if initially there is 
one unpaired atom inside the vortex at state \i). For this 
read-out scheme to succeed, the temperature needs to be 
kept as low as possible, so that there are no other un- 
paired atoms in the bulk around the vortex cores. How- 
ever, such thermally excited quasiparticles are expected 
to occur near the trap edges, and so for lasers sufflciently 
focused on the vortices near the center of the trap, one 
should be able to significantly suppress the unwanted 
signals. Finally, we mention that the resonant detec- 
tion laser may be replaced with a resonant multiphoton 
ionization process of the unpaired atom, yielding a sin- 
gle ionized electron that can be detected with essentially 
unit efficiency j^. Similar analysis as above also yields 
negligible impact on paired atoms. The advantage of the 
multiphoton ionization is that the detection may be done 
in parallel for all unpaired atoms in the sense that the 
electrons can be imaged on a channel plate, where the 
detection or no detection of the electrons is a parallel 
readout of all the qubits. However, such process is de- 
structive and the ionized atoms cannot be reused. 

In conclusion, we have proposed to use the two- 
dimensional, spin-polarized, Px + ipy atomic resonance- 
superfluid in the weak-pairing phase, potentially realiz- 
able in optical traps, in a suitable hardware for topo- 
logical quantum computation. We have given a realis- 
tic read-out scheme for the topological qubits, a major 
hurdle in this field, using the internal states of the con- 
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stituent atoms. 
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